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Skutterudites, a class of materials with cage-like crystal structure which have received considerable research 
interest in recent years, are the breeding ground of several unusual phenomena such as heavy fermion 
superconductivity, exciton-mediated superconducting state and Weyl fermions. Here, we predict a new 
topological insulator in bismuth-based skutterudites, in which the bands involved in the topological 
band-inversion process are d- and p-orbitals, which is distinctive with usual topological insulators, for 
instance in BiaScs and BiTel the bands involved in the topological band-inversion process are only 
p-orbitals. Due to the present of large d-electronic states, the electronic interaction in this topological 
insulator is much stronger than that in other conventional topological insulators. The stability of the new 
material is verified by binding energy calculation, phonon modes analysis, and the finite temperature 
molecular dynamics simulations. This new material can provide nearly zero-resistivity signal current for 
devices and is expected to be applied in spintronics devices. 

Topological insulator (TI) is a new kind of material which has gapped bulk state and gapless surface state with 
the latter protected by the topological character of TI' '°. For TIs with conserved spin along quantized axis, 
the topological order parameter is spin Chern number, and TI under time reversal symmetry is characterized 
by Zj quantum number''. The unique features of its surface state make TI have potential applications in spin- 
tronics and quantum information devices. TI is also the breeding ground for a good number of interesting 
quantum phenomena such as quantum anomalous Hall effect" Majorana fermions'^"' and topological mag- 
netoelectric effecf. TIs usually appear in those materials containing elements with strong spin-orbit coupling, for 
example, the bismuth element in Bi2Se3^'", BiTel"* '", and ScPtBi\ Moreover, pressure and strain has been 
demonstrated as an effective way to modulate the topological property of materials. For instance, CdSnAs2 under 
a 7% decrease in the lattice constant will become topological insulator^" while a 6% change in the length of c-axis 
will drive Bi2Se3 from topological non-trivial phase into topological trivial phase"". However, more interesting 
phenomena only can be induced by strong electronic interaction, such as the transition in correlated Dirac 
fermions^' and interaction induced topological Fermi liquids''''. Consequently, those TIs beyond p-band inversion 
arouse intensive research interest^^"^". Here, we predict a new d-p band-inversion topological insulator in bis- 
muth-based skutterudites in which the bands involved in the topological band-inversion process are d- and p- 
orbitals. Due to the present of large d-electronic states, the electronic interaction in this topological insulator is 
much stronger than that in other conventional topological insulators''^ 

Skutterudites, such as RhAss, IrAss, IrSbs and IrP3, crystallize in a cage-like crystal structure in which each 
transition metal atom octahedrally coordinates to six pnictide atoms"-''' (see Fig. 1 (a) where IrBi3 is illustrated). 
They have large Seebeck coefficients and therefore can behave as excellent thermoelectric materials''''. The 
discovery of heavy fermion superconductivity''^, exciton-mediated superconducting state'"' and Weyl fermions''^ 
in this system makes skutterudites a hot spot in condensed matter physics. Besides those skutterudites naturally 
exist, a number of new members in skutterudites have been experimentally synthesized, such as NiSbs*" in 2002 
and RuSb," in 2004. However, those materials are composed of elements with relatively weak spin-orbit coupling 
(SOC). Knowing that topological insulators are usually those materials containing elements with strong spin- 
orbit coupling strength, such as the bismuth element in topological insulator Bi2Se3^ '^, BiTel'"" and LaPtBi^'', It is 
reasonable to ask whether or not skutterudites composed of elements with strong spin-orbit coupling strength, i.e. 
bismuth, can exist stably and whether they can be topologically non-trivial? This new topological insulator in 
bismuth-based skutterudites, is exactly such kind of skutterudite material which is able to exist stably, contains 
elements with strong SOC, and has controllable topological phase transition. 

In this work, we predict a new d-p band inversion topological insulator in bismuth-based skutterudites, which 
is distinctive from usual topological insulators, for instance in Bi2Se3 and BiTel the bands involved in the 
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Figure 1 | Crystal structure and Brillouin zone, (a) unit cell of IrBis, including 8 Ir atoms (green) and 24 Bi atoms (pink) . Each Ir atom is surrounded by 6 
Bi atoms and each Bi atom has 2 Ir nearest neighbors, (b) The equivalent primitive cell of IrBi3, containing 4 Ir (green) atoms and 12 Bi atoms (pink), (c) 
The corresponding Brillouin zone and high symmetric points with T (0,0,0), H (0,1/2,0), N (1/4,1/4,0), P (1/4,1/4,1/4). (d) Free energy as a 
function of lattice constant. 



topological band-inversion process are only p-orbitals. Due to the 
present of large d-electronic states, the electronic interaction in this 
topological insulator is much stronger than that in other conven- 
tional topological insulators. The stability of the new material is 
verified by binding energy calculation, phonon modes analysis, 
and the finite temperature molecular dynamics (FTMD) simulations. 
We demonstrate that external strains are able to induce a topological 
phase transition in this system via band structure calculations. We 
confirm its topological non-trivial property by Zj quantum number 
calculation. 

Results 

Crystal structure and optimized lattice parameter. The bismuth- 
based skutterudite IrBis investigated here has space group IM3, and 
its crystal structure is shown in Fig. 1. There are 8 Ir atoms and 24 Bi 
atoms in a unit cell. Each Ir atom is surrounded by 6 Bi atoms and 
each Bi atom has 2 Ir nearest neighbors (see Fig. 1 (a)). The structure 
has space inversion symmetry with the inversion center (1/2,1/2,1/2). 
The structure belongs to the body-centered lattice type, and its 
primitive cell (Fig. 1 (b)) has a half volume of the unit cell. Fig. 1 
(c) shows the BriUouin zone and high symmetric points with T 
(0,0,0), H (0,1/2,0), N (1/4,1/4,0), P (1/4,1/4,1/4). 

We first optimize the lattice parameter and ionic positions. The 
calculated total free energy (solid line) as a function of lattice para- 
meter is shown in Fig. 1 (d). It can be clearly seen that the optimized 
lattice parameter (corresponding to the position of free energy min- 
imum) of the primitive cell is 8.493A. This value is 6% larger than 
that of IrSbs^^, which can be explained that Bi atom has a larger 
atomic radius than Sb atom. 



Binding energy calculation, phonon modes analysis and the finite 
temperature molecular dynamics simulations. In order to verify 
the stability of the new material, the authors perform the binding 
energy calculation, phonon modes analysis and the finite tempera- 
ture molecular dynamics (FTMD) simulations. The binding energy is 
calculated by 

Eb = EirBi, - nir -Eir - tlBi 'Ebi, ( l) 

where E/rSij denotes the free energy of IrBi3 per primitive cell, £/r and 
Egi the free energy of crystalline Ir and Bi per atom, rij^ and Wg, the 
number of Ir and Bi atoms in IrBi3 primitive cell. By simple 
calculation [There are «/r = 4 Ir atoms and n^i = 12 Bi atoms in 
an IrBis primitive cell. At GGA level, Ej^ = — 8.69 eV for crystalline Ir 
with space group FM3M and Egj = —3.70 eV for crystalline Bi with 
space group IM3M. From Fig. 1(d) we read Ei^Bi, — —82.81 eV. 
Substituting the above values in Eq.(l), we arrived at the binding 
energy = —3.65 eV.], Ei, is found to be equal to —3.65 eV per 
primitive cell. The negative value of binding energy infers a stable 
state of IrBis. 

Fig. 2 shows the phonon dispersion and phonon density of states 
(DOS) for IrBi3 at zero strain. In the phonon DOS subfigure, the 
black solid line represents the total phonon density of states, while 
the green and red shaded areas represent the states coming from Ir 
and Bi atoms, respectively. Phonon states in the low energy range are 
mostly composed of states from Bi atoms, indicating that Bi atoms in 
IrBi3 are much easier to vibrate than the Ir atoms. The phonon 
dispersion and phonon DOS show no imaginary frequency, indi- 
cating that IrBi3 is stable. 
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Figure 2 | Phonon dispersion and phonon density of states for liBi^. Orange dotted lines in all subfigures denotes the zero frequency. Calculations are 
performed at zero strain, (a) phonon dispersion curves for IrBis, in which the inset shows the dispersion near the zero energy, (b) phonon density of 
states for IrBij, in which black solid line represents the total phonon density of states, while the green and red shaded areas represent the states 
coming from Ir and Bi atoms, respectively. Phonon states in the low energy range are mostly composed of states of Bi atoms, indicating that Bi atoms in 
IrBis are much easier to vibrate than the Ir atoms. The phonon dispersion and phonon density of states shows no imaginary frequency, indicating that 
IrBis is stable. 



In addition, the dynamical stability of the material is further 
checked by finite temperature molecular dynamics simulations at 
temperature 300 K for room temperature and 30 K for low temper- 
ature. During the simulations, a 2 X 2 X 2 supercell containing 256 
atoms is used. The length of time-step is chosen as 5 fs and simula- 
tions with 1000 steps are executed. It is observed that, the atoms 
shake around the equilibrium positions back and forth while the 
extent of such motion under 300 K is larger than under 30 K (the 
evolution of atomic positions can be found in movies in supplement- 
ary information). However, no structural collapse happens through- 
out the simulations, which can also be seen from the free energies 
curves as the functions of time-step shown in Fig. 3. It is also 
observed that, the crystal structure always remains nearly the same 
as the initial crystal structure. Actually, as is shown in the inset of 
Fig. 3, the crystal structure corresponding to the last free energy 
maximum in T = 300 K case (right), still shows no significant struc- 
tural differences as compared with the initial crystal structure (left). 
The lattice relaxation, binding energy calculation, phonon modes 
analysis together with FTMD simulations mentioned above provide 
an authentic test for the stability of bismuth-based skutterudite IrBi3. 

Strain-induced d-p band-inversion topological insulator. The 

calculated band structures are listed in Fig. 4, where the black and 
blue lines represent the GGA and GGA+U band structures, 
respectively. As is shown in Fig. 4 (a), before exerting pressure, 
IrBia resides in the normal metal state with its bands crossing the 
Fermi level several times. Subfigure (b) to (d) represent the band 
structures at isotropic strain 3%, 6%, 9% respectively. With the 
increase of isotropic strain ((a) to (d)), the valence band crossing 
the Ep along H-N moves downwards and the density of states 



(DOS) at Fermi level decreases gradually. Under a 9% isotropic 
strain, the bands go across the Fermi level at F point but not at 
other points in BrUlouin Zone (BZ) (see Fig. 4 (d)), and the 
conduction band minimum and valence band maximum 
degenerate so that the material behaves as a semi-metal which 
have a zero energy gap, just like Graphene and CeOs4Asi2''. This 
degeneracy at F is protected by the cubic symmetry of crystal, 
which, as is tested by us, cannot be eliminated by small changes of 
the lattice constant. In order to shift the degeneracy at F, one needs to 
break that symmetry. An unsophisticated way is to add an anisotropy 
just like what was done on CdSnAs2^°. Here, we simply further 
impose a 2% suppression on the c-axis of the primitive cell while 
remaining the length of a- and b-axis unchanged, which imposes 
anisotropy on the system. While the anisotropy does not change 
the parities of each band, it opens a gap at the Fermi level, 
dragging the system in the insulating state (see Fig. 4 (e)). Fig. 4 (f) 
shows the Ir-d projected band structure near the Fermi level and near 
F point, in which the radii of red circles correspond to the proportion 
of Ir-d electrons. It can be seen that, those localized bands above the 
Fermi level are mainly contributed by d-orbitals of Ir atoms. The 
highly dispersive band below the Fermi level is mainly contributed by 
p-orbitals of Bi atoms, and it has little weight of Ir atoms in those k- 
points far away from F point. However, in the vicinity of F point, the 
weight of Ir atoms in that band increases rapidly and becomes 
dominating orbital component, showing an apparent band 
inversion. Such band-inversion character is further checked by the 
modified Becke-Johnson (mBJ) potential (see supplementary 
information), which is proved to be able to predict an accurate 
band gap and band order'"''"^''. In order to further confirm the topolo- 
gical property in such condition, we calculate the Zj topological 
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Figure 3 | Finite temperature molecular dynamics. Free energies as functions of time-step at temperature T = 30 K (blue curve) and T = 300 K (red 
curve). The slight shift of the free energy curves corresponds to the oscillations of each atom around their equilibrium position. The absence of sharp 
changes in such curves indicates that no structural phase-transition happens throughout the whole simulation process. The initial crystal structure 
(denoted by the orange circle on the free energy curve) is plotted in inset (a). The crystal structure corresponding to the last free energy maximum 
(denoted by the green circle on the free energy curve) is shown in inset (b) as a comparison. It can be seen that, the latter still shows no significant structural 
differences as compared with the initial crystal structure. 
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Figure 4 | Band structures of IrBis. The black and blue lines in all subfigures represent the GGA band structures and GGA+U band structures 
respectively, (a) band structure without exerting pressure, the system is in normal metal state with its bands go across the Fermi level several times, (b) to 
(d) represent the band structures at isotropic strain 3%, 6%, 9% respectively. With the increase of isotropic strain ((a) to (d)), the valence band crossing 
the £p along H-N moves downwards gradually. In the band structure under 9% uniform strain (d), a zero gap metal state is obtained, (e) further impose a 
2% suppression on the length of c-axis of the primitive cell, a gap appeared at the Fermi level due to the breaking of the cubic symmetry. The inset of (e) is 
the zoom-in of the band structure close to the Fermi level, (f) Ir-d projected band structure near Fermi level, the radii of red circles are proportional to the 
weight of Ir-d states, showing a significant band inversion. 
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Table 1 | Paritiesof top-most isolated valence bands at eight time-reversal invariant momenta. Positive parity is denoted by -I- while negative 
denoted by — . Products of the occupied bands at each time-reversal invariant momentum are listed in the right-most column. As is shown, the 
productof paritiesof occupied bands contributes a —1 at (0,0,0) while + 1 at the seven other time-reversal invariant momenta, resulting in vo 

= 1, V, = V2 = V3 = 0 



(0,0,0) + + + + + - + + 

(7t,0,0) - + - + - + - + - + + + + 

[0,11,0) + - + + + + + - + - + - + + + + 

[0,0,n) - + - + + + + + + 

(7t,7t,0) + - + + + + + 

{ti,0,ti) - + - + - + - + - + + + + 

[0,11,71) + - + + + + 

[tz,i:,k) + + + + + + + - + + + + + - + 



quantum number of the system by the Fu-Kane method''. The index 
for strong topological insulators Vg is expressed as ( — 1)''° = Yl^^ j 5,- 
in which i5,= = i iimi^i) represents the product of the parities of 
the occupied band at 8 time-reversal invariant momenta T,. The 
calculated parities of top-most isolated valence bands (here refers 
to the isolated block of states between —8.0 eV to 0 eV in Fig. 5) 
at eight time-reversal invariant momenta are listed in Table 1, where 
the deeper states (those states lower than —9.5 eVinFig. 5) separated 
far from top-most isolated valence bands are ignored because they 
don't change system's band topology. As is shown, the product of 
parities of occupied bands contributes a — 1 at F while + 1 at the 
seven other time-reversal invariant momenta. As a result, Z2 
quantum number is Vg = 1, Vj = V2 = V3 = 0, which corresponds 
to a strong topological insulator. 

Partial-density of states and the d-p orbitals dominating property 
near the Fermi level. Fig. 5 depicts the atomic- and orbital-resolved 
density of states (DOS). The black solid lines in subfigure (a)(b)(c) 
represent the total DOS. Fig. 5 (a) is atomic- resolved DOS, in which 
the green curve represents the states of Ir while the red curve 
represents the states of Bi. It's clear that the DOS of both types of 
atoms is in quite large values, indicating that both types of atoms 
make a significant contribution to the total DOS. This is different 
from M0S2 where states near Fermi level are dominated by only one 
kind of atom (Mo)^''. Fig. 5 (b) and (c) are orbital- resolved DOS of Ir 
and Bi atoms, respectively. Green, blue and red curves represent s-, p- 
and d-orbitals. One character of the material introduced here is a 
large proportion of d-states near Ep. 

Further, we calculate the d-orbitals projected band structures (see 
Fig. 6) in the local coordinate of the Bi octahedral. The orange, violet, 
red, green and blue colors in Fig. 6 represent the d^i , (1^2 ^yi , d^y, dy^ 
and dxz orbitals respectively. The radii of circles are proportional to 
the weights of corresponding orbitals. It can be seen that, the 
orbitals (including the d^y, dy^ and d^^ orbitals) reside far below the 
Fermi level and are fully occupied. While, the lowest three conduc- 
tion bands are mainly contributed by the e^ orbitals (including the d^2 
and dx^^yi orbitals). More specifically, the dj.2__^2 orbital makes an 
even larger contribution than the the d^i orbital for the lowest con- 
duction band. The large proportion of d-states near Ep is distinctively 
different from usual TI materials, for example states near Ep mainly 
containing s-p electrons in HgTe and p-electrons in Bi2Se3. The large 
proportion of d-states near Ep indicates that electrons in such mater- 
ial process strong electronic correlations and are more localized than 
other TI materials. The strong electronic correlations make the 
material a good platform for investigating the effect of correlations 
on the topology, as well as a candidate for realizing the quantum 
information device based on correlations. The localization will 
enhance the effective mass of bulk electrons, and hence reduce the 
bulk contribution to the local current at finite temperature, making 
the spin-binding property more apparent, which is helpful in fabric- 
ating spintronics devices with higher stability. 



Discussion 

Experimentally, the new strain- induced topological insulator IrBi, 
could be grown using the Bridgman method, by which the C0P3" 
and the RuSba^* crystals have been successfully synthesized. The 




Energy (eV) 

Figure 5 | The atomic- and orbital-resolved density of states. The black 
solid lines in all subfigures represent the total density of states (DOS), (a) 
atomic-resolved DOS, in which the green curve represents the states of Ir 
and the red curve represents the states of Bi. It's clear that both type of atom 
made a significant contribution to the total DOS, different from M0S2 
where states near Fermi level are dominated by Mo. (b) and (c) are orbital- 
resolved DOS of Ir and Bi atom respectively. Green, blue and red curves 
represent s-, p- and d-orbitals. 
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Figure 6 | Orbital-projected band structures. The orange, violet, red, green and blue colors in subfigures represent the d^i , d^i _yi and A^y, dy^ and 
orbitals respectively. The radii of circles are proportional to the weights of corresponding orbitals. The Fermi level is set to be zero energy. It can be seen 
that, the orbitals (including the d^y, dy^ and d^^ orbitals) reside far below the Fermi level and are fully occupied. While, the lowest three conduction 
bands are mainly contributed by the e^ orbitals (including the d^i and d^^-f orbitals). More specifically, the dj^i^yi orbital makes an even larger 
contribution than the the d^i orbital for the lowest conduction band. 



crystal growth should be conducted in a sealed quartz ampoule. The 
iridium and bismuth should be coated by graphite and then intro- 
duced into the quartz ampoule. A temperature gradient of about 
50°C/cm should be maintained at the growth interface, just like in 
the case of RhSb3'^. To remove the excess bismuth in the as-grown 
crystal, post-annealing should be performed'''. After the synthesis of 
the new material, its crystal structure could be characterized by the X- 
ray diffraction using the monochromatic Cu Ka radiation'^. Then, 
the strains could be generated by a pair of diamond anvils"*, which 
was used to generate strong pressure even above 200 GPa''. 
Moreover, the real-time pressure strength could be detected by ruby 
fluorescence method"*'"'. In order to verify the topological property 
of the material, it is suggested to perform the transport measure- 
ments'". Similar to 612863, the observation of the spin-Hall currenf*^ 
and the non-equally spaced Landau levels'*'' in IrBis will be signatures 
of the Dirac fermions in surface of the topological insulator'. 

In this work, we predict a d-p band inversion topological insulator 
bismuth-based skutterudite IrBi3, and verify its stability. Our results 
indicate that this material is zero gap semi-metal after imposing 
uniform strain, and it can become topological insulator if an aniso- 
tropy is further applied to break the cubic symmetry. Furthermore, 
near the Fermi level there is a large proportion of d-electronic states 
which is distinctive from usual topological insulators, for instance in 
612863 and BlTel the bands Involved In the topological band-inver- 
sion process are only p-orbitals. Consequently, the electronic inter- 
action in this topological insulator is much stronger than that in 
other conventional topological insulators. This provides realistic 
material for investigating the effect of correlations on the topology, 
fabricating quantum information devices and spintronics devices 
with higher stability. 

Methods 

Our first principle calculations are in the framework of the generalized gradient 
approximation (GGA) of the density functional theory. The VASP package'^^-^" has 
been employed and the projector-augmented-wave pseudo-potentials^^ are used. 
Plane waves with a kinetic energy cut-off of 400 eV are used as basis sets and k- 
point grids in BrUlouin zone is chosen as 6 X 6 X 6 according to the Monkhorst-Pack 
scheme. The relaxations are carefully made so that the forces on atoms are smaller 
than 0.0003 eV/A, in which the conjugate gradient algorithm is utilised. In the finite 



temperature molecular dynamics simulations, a 2 X 2 X 2 supercell containing 256 
atoms is used and the length of time-step is chosen as 5 fs. The phonon dispersion 
curve and phonon density of states are obtained using the force-constant method by 
phonopy code^^. The effect of spin-orbit coupling (SOC) is included in the calcula- 
tions after the structural relaxations. GGA-t-U calculations are based on the 
Dudarev's approach implemented in VASP, with the effective on site Coulomb 
interaction parameter U — 3.0 eV and the effective on site exchange interaction 
parameter / — 0.5 eV for d-orbitals of Ir atoms^**. The GGA band structures are 
checked by the full-potential DFT code WIEN2k^^ in the supplementary information. 
We also use the modified Becke- Johnson (mBJ) semilocal exchange-correlation 
potentiaP'''^'' to further check the band order and the magnitude of energy gap, and in 
this process the GGA wave function is used to initialize the mBI calculation. 
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